Brownian Entanglement. 



O 

o 



A.E. Allahverdyan^'^\ A. Khrennikov^' and Th.M. Nieuwenhuizen^^ 
Institute for Theoretical Physics, Valckenierstraat 65, 1018 XE Amsterdam, The Netherlands 
Yerevan Physics Institute, Alikhanian Brothers Street 2, Yerevan 375036, Armenia and 
^^International Centre for Mathematical Modeling in Physics and Cognitive Sciences, 
MSI, University of Vdxjd, S-35195, Sweden 
(Dated: Today: February 1, 2008) 

We show that for two classical brownian particles there exists an analog of continuous-variable 
quantum entanglement: The common probability distribution of the two coordinates and the cor- 
responding coarse-grained velocities cannot be prepared via mixing of any factorized distributions 
referring to the two particles in separate. This is possible for particles which interacted in the past, 
but do not interact in the present. Three factors are crucial for the effect: 1) separation of time- 
JnI ^ scales of coordinate and momentum which motivates the definition of coarse-grained velocities; 2) 

the resulting uncertainty relations between the coordinate of the brownian particle and the change of 
its coarse-grained velocity; 3) the fact that the coarse-grained velocity, though pertaining to a single 
brownian particle, is defined on a common context of two particles. The brownian entanglement 
O ■ is a consequence of a coarse-grained description and disappears for a finer resolution of the brow- 

CD ■ nian motion. We discuss possibilities of its experimental realizations in examples of macroscopic 

I brownian motion. 

^ ! PACS numbers: 05.40.Jc, 03.65.Ta, 05.70.Ln 

^ . I. INTRODUCTION. 

(N 

CO ' There is a long tradition of seeking connections between quantum mechanics and classical statistical physics. A list 
of known examples includes i) analogies between the notion of complementarity in quantum mechanics and statist!- 
cal thermodynamics, in particular, between quantum mechanical uncertainty relations and the energy-temperature 
uncertainty relation in statistical thermodynamics ii) mathematical relations between Schrodinger and 

Fokkcr-Planck equations 0, which makes quantum intuition very useful, e.g., for polymers 0; in) attempts to derive 
Schrodinger equations from the classical kinetic picture or from classical stochastic electrodynamics 

More generally, both quantum mechanics and statistical physics are essentially probabilistic theories and already at 
this level one expects to find certain analogies between their concepts. The progress in this direction has been for a long 
^ _ time plagued by statements on incompatibility between quantum mechanics and classical probability theory, which 
^ • lies in the basis of classical statistical physics. As an example see Ref. iQj for statements that quantum interference 
^ ' phenomena contradict to classical probability theory. It was, however, noted that such opinions are not warranted 
[TOj, and recently it was shown explicitly that several basic relations of quantum mechanics can be derived from 
the classical probability theory provided the contexts of physical conditions (measurements) are properly taken into 
account [TH . 

Once the basic probabilistic ground of quantum mechanics and classical statistical physics is recognized to be the 
] same, one wonders whether some unusual aspects of quantum mechanics, such as entanglement, can find analogies in 
classical statistical physics. We see at least two reasons for seeking such analogies. First, it is going to refresh our 
understanding of classical statistical physics, and may imply in future that advantages offered by quantum mechanics 
in certain tasks of information processing and transfer are not unique to quantum mechanics, and can be looked 
for in classical areas of physics as well Second, it is useful for understanding quantum mechanics that certain 

concepts believed to be purely quantum — that is, incomprehensible in classical terms — can find natural classical 
analogies. In fact the works and ideas mentioned in the above point were partially directed toward this goal. 

More recent results along these lines are classical analogies to quantum entanglement found in optics |^ and 
classical information theory ,23i| (secret classical correlations), and a classical probabilistic model for certain aspects 
of (multi-time) quantum measurements [23| . 

The purpose of the present paper is to show that quantum entanglement can have a natural analogy in the physics 
of brownian particles (brownian entanglement). The reason for the existence of this analogy can be qualitatively 
explained as follows. It is known that the dynamics of a brownian particle can be observed at two levels 5]. Within 
the first, more fundamental level the brownian particle coupled to a thermal bath at temperature T is described via 
definite coordinate x and momentum p and moves under influence of external potential, friction force and an external 
random force. The latter two forces are generated by the bath. The second, overdamped regime applies when the 
characteristic relaxation time of the coordinate is much larger than that of the momentum Tp, Tp (overdamped 

regime). On times much larger than Tp one is interested in the change of the coordinate and defines the coarse-grained 
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velocity as u = Ax/ At for ^ At ^ Tp. This definition of v is the only operationally meaningful one for the 
(effective) velocity within the overdamped regime. It appears that the coarse-grained velocity, though pertaining to 
single particles, is defined in the context of the whole systems of coupled brownian particles. Together with uncertainty 
relations between the coordinate and the change of the coarse-grained velocity — the role of Planck's constant is being 
played by the temperature of the bath, — this contextuality feature will be shown to cause a phenomenon similar to 
quantum entanglement: The common probability distribution of the two coordinates and the corresponding coarse- 
grained velocities — for two brownian particles which interacted in the past, but need not interact in the present 

cannot be prepared via mixing of any factorized distributions referring to the two particles separately. This brownian 
entanglement is a consequence of a coarse-grained description and disappears within the first (more fundamental) 
level of description, simply because entanglement is absent in classical mechanics. 

The paper is organized as follows. In section ^ we recall the phenomenon of (continuous- variable) quantum 
entanglement, focusing especially on the relations between the entanglement and the uncertainty relations. Section 
mil discusses coarse-grained velocities and uncertainty relations for the classical brownian motion. Next two sections 
define and study the phenomenon of brownian entanglement. Sections Ivl and IvTl offers detailed comparison between 
the features of quantum entanglement and those of its brownian counterpart. In section lVlIl we discuss possibilities of 
experimental realization of brownian entanglement. Our conclusions are presented in the last section. Some technical 
questions are worked out in Appendix. 

II. QUANTUM ENTANGLEMENT 
A. Statistical interpretation. 

This section recalls the phenomenon of entanglement in quantum mechanics, and especially underlines its connec- 
tions with the uncertainty relations. 

Before starting, it is useful to stress that in the present paper we adhere to the statistical (ensemble) interpretation 
of quantum mechanics, where a quantum 'state' is described by a density matrix p, and any state, inclu ding a pure 
state refers to an ensemble £{p) of identically prepared systems; see, e.g., [l^j llJj llM llfl IISj llli IS ^- To 

put it succinctly: quantum mechanics makes^ 

• statistical statements on 

• the results of measurements done 

• on ensembles of identically prepared systems. 

As was stressed repeatedly 0, 0, 0, 0, 0, 0, , in particular by experimentalists , the experimentally 
relevant statements of quantum mechanics do not require more than the minimal statistical interpretation Moreover, 
this interpretation deals more successfully with the measurement problem, as instanced by a recent exactly solvable 
model '27] , and allows to reconcile quantum mechanics with classical probability theory . The fact that discussions 
on quantum entanglement (and on Bell inequalities and related matters) do not employ the statistical interpretation 
is a mere prejudice; see Refs. 0,|25, 26, 27] for examples of such discussions. 

B. Definition of entanglement. 

Consider a quantum system S consisting of two subsystems Si and S2. A state p of S is called entangled, see e.g. 
|12|. (or non-separable) with respect to Si and S2, if it cannot be represented as 

n n 

fc=i fc=i 



^ The minimal statistical interpretation should of course be distinguished from various hidden- variable theories and assumptions, in 
particular, from the pre-assigned initial value assumptions, where quantum measurements are viewed as merely revealing pre-existing 
values of all observables. Unfortunately, some proponents of the statistical interpretation were unc lear at this point, a fact that for a 
while discredited this interpretation. For a well-balanced discussions on this and related points, see |l8lllflll . 

^ This way of putting the message of stati stical interpretation is adopted from l26i . 

^ Still, as correctly pointed out in |2,'ill3l| . the choice of interpretation can and does influence one's estimates of importance for various 
scientific problems. 
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where p^^^ and pf'^ are arbitrary density matrices living in the Hilbert spaces of Si and S2, respectively, n is an 
integer, and where {pk\^=i is a probability distribution. 

According to definition J^lj a separable quantum state can always be prepared by means of mixing ^ non-correlated 
states /5^^^ (g) p^^^ of the two subsystems. For a pure state p — we return to the more known definition of 

entanglement: IV') cannot be represented as = |'0)i ® IV')2- For this particular case the absence of entanglement 
implies the absence of any correlation. In contrast, for the more general situation given by a non-entangled state 
can still posses certain (classical) correlations, since the totally uncorrelated situations will be given as p = p'^^ ® p'^'^\ 

Two important features of quantum entanglement should be noted. First, it can be tested only through measuring 
some correlations between the subsystems Si and S2. Observables pertaining to Si or to S2 alone are obviously 
insensitive to entanglement. Second, if the state is prepared by two different observers 1 and 2, then this preparation 
is seen to involve correlated actions of them, and thus the corresponding observers have to communicate classically. 

C. Classical systems. 

In classics the representation is apparently always possible. Indeed, let us have a probability distribution 
P{xi,X2) of two classical systems Si and S2 represented by random variables xi and X2- Assume for simplicity that 
Xi can take values Oi, -.jan, while X2 can take values bi, .., 6„. Then P{xi, X2) can be written as 

Pix,,x2)= E p(^^0)^ 

axi ^X2l3: (2) 

Q=ai,...,a„ /3=6i,...,6„ 

where Saf3 — Ufa = (3, and Safj ~ otherwise. Provided that 5axi and 6x2/3 are two legitimate probability distributions 
belonging to Si and S2, that is, provided there are no mechanisms prohibiting the realization of 5axi and Sx2i3 as 
physically acceptable distributions usable in the actual preparation, the representation ^ is realized for the classical 
situation: there is no entanglement. 

In fact the quantum situation goes to the classical one right in if we assume that the involved density matrices are 
always diagonal. So it is the presence of non-diagonal elements of the given density matrix that makes the situations 
different. 



D. Semiclassical systems. 

We shall now localize the cause of quantum entanglement for systems which behave (semi) classically in several 
other respects. This will help us to understand the way of searching for entanglement in non-quantum situations. 

Let Si and S2 be two quantum (non- interacting) harmonic oscillators, and the overall density matrix (in the 
coordinate representation) p{xi, x'^; X2,x'2) be a gaussian function of its variables. The corresponding coordinate and 
momenta operators are denoted by xi, pi and £2, P2, respectively. As well-known, gaussian states are conveniently 
dealt with help of Wigner functions W{x,p), which are functions of x and p ^. The Wigner function is equivalent 
to the density matrix, that is, it allows to calculate all averages, but in contrast to the density matrix it has several 
properties expected for the common probability distribution of the coordinate and momentum. In particular, the 
analysis via Wigner functions will provide us below with a richer intuition on the relations between entanglement and 
uncertainty relations; thus, between entanglement and non-commutativity. 

Since the overall Hamiltonian of Si and S2 is assumed to be harmonic, the initially gaussian state remains gaussian 
for all times and the corresponding Wigner function will be positive and thus will partially admit a classical interpre- 
tation in terms of common probability distribution of the coordinate and momentum '3j|. A related fact is that for 
two harmonic oscillators the Ehrcnfcst equations of motion can be recast into a classical form |i33i] . So many aspects 
of this system can be accounted for in classical terms. 



* Mixing ensembles £{pi) and £(p2) with probabilities pi and p2, respectively, means that one throws a dice with probabilities of outcomes 
equal to pi and P2, and depending on the outcome one picks up a system from £{p\) or £{p2), keeping no information on where the 
system came from. Alternatively, one can join together Npi systems from £{pi) and Np2 systems from £{p2) {N ^ 1), so that no 
information information is kept on where a single system came from. 

^ For good reviews on the properties of Wigner function see l33l . 
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In order to see why in spite of these classical features the system of two oscillators can be entangled, note that for 
the Wigner functions one can rewrite the condition ^ as 

Wixi,pv,X2.P2)= J dXr{X)Wi{xi,pi\X)W2ix2,P2\X), V{X)>0, JdXV{X) = l. (3) 

where Wi and W2 are separate Wigner functions for Si and S2 respectively, and V{X) is some probability distribution. 
If Si and S2 were classical oscillators, then instead of Wigner functions we would have distribution functions, and one 
can always write down the analog of 10) for the common distribution function P(xi,pi; X2,P2) and A = (ai 012, (32)'- 

P{xi,pi;x2,P2) = y"'i"i "^/^i da2d(32 P{ai, Pi; a2, (32) Ho^i - xi)S{(3i - pi) 6{a2 ~ X2) 5{P2 ~P2)- (4) 

In the classical situation this means that there is no entanglement. The same formula (0J can formally be written 
down also for the positive Wigner functions W{xi,pi]X2,P2)- However, in the quantum situation 5{ai —xi) S{Pi —pi) 
and d{a2 — X2) S{(32 — P2) are not legitimate Wigner functions, since they prescribe definite values to both coordinate 
and momentum and thus do not respect the uncertainty relations. 

We conclude that the uncertainty relations are necessary — but not sufficient — for the existence of entanglement 
in semiclassical systems. 



E. A simple sufficient condition for quantum entanglement 

The definition of the entanglement as given by Eq. ^ is not practical (except for a pure density matrix /5, when 
no entanglement means no correlations). It is therefore useful to have certain sufficient conditions for the presence of 
entanglement which will be easy to handle in applications and which will have a transparent physical meaning. 

We choose the units in such a way that the coordinate and the momentum have the same dimension as ^/h. For 
a harmonic oscillator with mass m and frequency w it will suffice to make the following canonical transformation: 
X — !■ ^muj X, p — > pI ^vtill}. 

For a two-particle system S with coordinate and momenta operators xi, X2 and pi, p2, respectively, one can propose 
the following sufficient condition for the entanglement p^ . 
Let us first note that the standard uncertainty relation 

(Ai^) (Ap^) > ^, (5) 

Ai; = i;— (f), Ap = p— (p), (6) 

implies 

(Ai2) + (Ap2)>(Ai^) + ^^, (7) 

and then one gets via minimizing the RHS of {Tj) over (Aa;^) (which produces (Ai^) ft/2, to be put in the RHS of 
0): 

(Ax^) + (Ap2) > h. (8) 
Now assume that the two-particle system is described by a factorized density matrix: 

p = (9) 

Then due to lO one can write: 

((Ail - Ax2)') + ((Api + Ap2)') = (A£?) + (Ap?) + (Ax^) + (Ap^) > 2ft, (10) 

just because for non-correlated systems the corresponding variances add up. If the overall system is in a separable 
state, then ifTUIl is even strengthened, since the variance of any observable (AA^) increases under mixing, i.e. under 
a transformation: 

{Pk-,Pk\ ^^Y^Pkh-, Pk>0, ^pfc = 1, (11) 
k k 
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where pk are normahzed (trp^ = 1) density matrices. Indeed, 

(Ai') = ^ Pfctr [pk{A-{A)f] = ^ pfe {tr [ pfc 12 ] - tr [ pfc i } + ^ Pk {tr [ p,. i ] - (i) }' 

= ^ Pfc (Ai^) + ^ { tr [ pfe i ] - (i> } ' > ^ p,. ( Ai^),.. (12) 

Employing p2|) for A = Axi — Ax2 and for A = Api — Ap2, we get that though (|10l) was obtained for the factorized 
state (O, it remains valid for an arbitrary non-entangled state J^kPkp'k^ ® p'k'' ■ Thus the violation 

{{Axi - Ax2f) + ((Api + Ap2f) < 2h (13) 

of p(J|) is a sufficient condition for entanglement. 
Eq. (|13|l has a transparent physical meaning: 

• entanglement is present, if fluctuations do not sum up additively, i.e., if the changes Axi and Ax2 of the 
coordinates tend to correlate with each other, while those of the momenta tend to anticorrelate. 

Alternative sufficient conditions for entanglement can be built based on {{Axi + Ax2)'^) + {{Api — Ap2Y) or 
((Ail + AX2Y) + ((^Pi + Ap2)2), etc. Various conditions obtained in this way are obviously not equivalent to each 
other. An obvious way to strengthen it is to demand that 

{{Axi + eAx2f) + ((Api + CAp2)') < 2ft, (14) 

at least for one of four independent choices e = ±l, C = ±l- 

For some special states (e.g. gaussian states) there exist in literature necessary and sufficient conditions for en- 
tanglement [30j. However, they are technically involved and do not have such a straightforward physical meaning as 
(I13|l . For our purposes conditions (|13l I14|l are sufficient. 



F. Operational issues. 

Let us finally recall how the condition (|13|l is checked operationally. To this end, we re-write this inequality as 

{AxD + {AxD + {ApD + {ApD (15) 
-2(Ai;i A2;2) +2(Api Ap2) (16) 

< 2h. 

One now needs two ensembles of systems S, each one consisting of identically prepared correlated sub-systems Si 
and S2. Assuming that Si and S2 are in possession of observers 1 and 2, respectively, the observer 1 measures on 
the first (second) ensemble xi (pi), while the observer 2 measures on the first (second) ensemble X2 (^2)- Note that 
the above two ensembles are necessary, since xi and pi (respectively, X2 and P2) do not commute: [xk,pi\ — ifiSki, 
with Ski being Kronecker's delta, and thus cannot be measured simultaneously. Now the quantities in 115|) can be 
estimated by each observer separately, while for the quantities in IjlGI) the observers need to put the results of their 
measurements together (or to communicate in any other classical way) and to count the coinciding events. 



III. COARSE-GRAINED VELOCITIES AND UNCERTAINTY RELATIONS FOR BROWNIAN 

PARTICLES. 



Consider N identical brownian particles with coordinates x = (xi, xjv) and mass m interacting with N indepen- 
dent thermal baths at temperatures Ti and subjected to a potential U{xi^ ...,xn). The overdamped limit is defined 
by the following two conditions '3| (for a more detailed discussion see Appendix): 

• The characteristic relaxation time ^ of the (real) momenta mi, is much smaller than the one of the coordinates. 



The characteristic relaxation time Tg of a random time-dependent variable 6 (e.g. coordinate or momentum) is an ensemble notion and 
is defined as the time necessary for the conditional distribution of to become memoryless: P{d,t' + r0\6' ,t') = P{0,t' + Tg). For the 
system of brownian particles sufficiently strongly coupled to their thermal baths, the arbitrarinesses of the above definition (e.g., the 
precise choice of 9' and t') do not change the physical content of the definition . 
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• One is interested in times which are much larger than the relaxation time of the momenta, but which can be 
much smaller than or comparable to the relaxation time of the coordinates. 

Under these conditions the dynamics of the system is described by the following Langevin equations |^: 

= /,(x)+ry,(i), /,(x) = -a,.(7(x), {r^,{t) rjj{t')) ^ 2T, 6,, S{t - 1'), (17) 

where for our convenience both the mass m and the damping constant (coupling constant of the particles to the bath) 
are in the main text taken equal to one (they are recorded in the Appendix) ''. 

The conditional probability /'(x, i|x', t') is known to satisfy the following Fokker-Planck equation 

dtP{^, i|x', i') = - E W ^(^' ^1^'' ] + E ^1^'' ^'), t > t'. (18) 

i i 

Eq. l(T^ is associated with the following initial condition: 

JV 

P(x, t|x', t) = (5(x - x') = n S{x, - x',), (19) 



as follows from the very definition of the conditional probability. 

Consider an ensemble I](x, t) of all realizations of the whole A^-particle system which at time t have a coordinate 
vector X. Such an ensemble can be selected out of all possible realizations by measuring x — {xi, ...,xn)- For this 
ensemble the average coarse-grained velocity for the brownian particle with index j might naively be defined as: 

«,(x,t) =lime^o J dy ^'"^' P{y,t + e\^,t). (20) 

However, the absence of regular trajectories enforces us to define different velocities for different directions of time 01 : 

w+,j(x, t) = lim,_,+o J dy, yi—^ P{y^^t + e|x, t), (21) 

z;_,,(x, t) - lim,^+o J dy, ^1—^1 P{y^,t ~ e|x, t). (22) 

The physical meaning of these expressions can be explained as follows. 

1. As seen from the definitions, w+.j(x, t) is the average velocity to move anywhere starting from (x, i), whereas 
?;_.j(x, t) is the average velocity to come from anywhere and to arrive at x at the moment t. 

Since these velocities are defined already in the overdamped limit, e is assumed to be much larger than the 
characteristic relaxation time of the (real) momentum which is small in the overdamped limit. Therefore, 
we call (|21l I22II coarse-grained velocities. It is known that for the overdamped brownian motion almost all 
trajectories are not smooth. This is connected to the chaotic influences of the bath(s) which randomize the 
real momenta on much smaller times, and this is also the reason for u_|_.j(x, t) ^ w_j(x, t). The difference 
w_i- j (x, t) ~ u_ j(x, t) thus characterizes the degree of the above non-smoothness. 

Recall that would one take e much smaller than the characteristic relaxation time of the momentum — which 
will amount to applying definitions (|21(l and H22|) to a smoother trajectory — then w+j (x, i) and v_.j(x, i) will 
be equal to each other and equal to the average momentum. These points are discussed in detail in Appendix. 



Eq. 1171 can be obtained from the complete Langevin equations of motion: [with the same definition of 'r]i(t) as in I17i ] mxi + Xi = 
/i(x) + Tjiit), by disregarding the first term mXi corresponding to acceleration. One might make this heuristically by requiring that the 
friction force Xi dominates, and that the second-derivative Xi is small for long times. More rigorous derivation of this classical problem is 
presented in 0, and is recalled in Appendix for a simple model. In their turn the above complete Langevin equations can be rigorously 
derived from the Newton equations of motion for the brownian particles and the baths. Here are the basic conditions for this derivation. 
1) The thermodynamical (macroscopic) limit for the baths. 2) A reasonable model for the particle-bath interaction. 3) The initially 
equilibrium state of the baths, described by the corresponding Gibbs distribution s at te mperature Ti (this is how the random noises 
rii(t) come into existence). The details of the derivation can be looked up, e.g., in mil 
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2. Here is an operational procedure for measuring v+j{'x,t) and w_j (x, i). Consider the overall ensemble of the 
brownian particles. The single members of this ensemble are N coupled brownian particles. For each such single 
member one measures 

i) the coordinate of the brownian particle with index j at the moment t — e; 

ii) the coordinates of all brownian particles at the moment t; 

Hi) the coordinate of the same brownian particle with index j at the moment t + e. 

Repeating these points many times on various members of the above overall ensemble, ignoring all the results 
from i) and in), and selecting events from the second step, we reconstruct the ensemble S(x, t). Employing 
the results from i) and Hi) and conditioning repeatedly upon a single member from E(x, <), we estimate the 
conditional probabilities P(yj, < ± e|x, t) and then finally calculate w+j(x, and u_j(x, via Eqs. H21II22|I . 

3. The above ensemble S](x, t) is different from an ensemble Yij{xj,t) obtained by measuring Xj (at time t) irre- 
spective of other coordinates. Indeed, S(x, t) is defined globally, and if there are different observers in possession 
of each brownian particle, they have to communicate to each other in order to be able to construct I](x, <). In 
contrast, Y,j{xj,t) is defined exclusively with respect to the brownian particle with index j. This point will be 
discussed in more details later on in section Ivl 

The calculation of v+^i{x,t) and v-^i{x,t) is straightforward upon using the following three things: relation 

F(x, t + e|y, t)=S{^-y)+eJ2 [-My)d., ^(x - y) + T, dl,^5{^ - y)] , (23) 

i 

which follows from 1)181 119|l , Bayes formula, and partial integration assuming natural boundary conditions at infinity: 
P{x, i) ^ if Xj — > ±oo. Starting from the definitions H21II22|) . we obtain 

v+,j{x,t) = lim,^+o J dy^^^P(y,i + £|x,<) 

= 1 dy {y, -x,)Y,[ -ni^)dyA^ - y) + T..dlyA^ - y) ] 

i 

= /j(x), (24) 



z;_j(x,i) = lim,^+o J dy^^-^P{y,t~e\^,t) 
= Imie^+o / dy ' ^ P(x,%,t - e) 



P(x,i) 

^""^ 1^) ^ ["^^^^ ^^^-'^^'^ " + T,d,..A^ - y)] 

= /,(x)-2r,9,JnP(x,i). (25) 
The difference between the coarse-grained velocities v+j{x,t) and u_j (x, i) is 

^^.(x,i) = ^-A^,t)~^+A^^t) ^ InPM. (26) 

Recall that this quantity is non-zero due to the action of the thermal bath (the factor Tj in p6() ). and due to the 
fact that at the coarse-grained level of description almost all trajectories of the brownian particles are not smooth. 
Sometimes uj (x, t) is referred to as osmotic velocity of the brownian particle with index j Q . We shall use this word 
as a useful short-hand for the more precise term "change of the coarse-grained velocity" . 

Note that once the interactions between the brownian particles are absent: U{x.) = U{xk) — and consequently 
correlations are absent as well, P(x, t) = P{xk, t), if they were absent initially — the coarse-grained velocities w± j 
depend only the corresponding coordinate: z;±j(x, i) — v±j{xj,t). 

The regular (or mechanical) counterpart of the coarse-grained velocity for the brownian particle with unit mass 
and index j can be naturally associated with the newtonian force /j(x) as defined by H17() . When the bath(s) are 
absent, this is the only contribution to the coarse-grained velocity. The fact that velocity appears to be proportional 
to force reminds us that we are in the overdamped regime of description, where, in particular, friction dominates over 
acceleration. 
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A. Uncertainty relations and their interpretation. 

Uncertainty relations exist not only in quantum mechanics but also in physics of overdamped brownian motion |29j| . 
As the very subject of statistical physics, they arise out of ignorance reasons, or more precisely due to separation of 
time-scales: though the real classical particles involved in the classical brownian motion certainly do have sharply 
defined coordinates and momenta at the microscopical level, at the coarse-grained (overdamped) level of description 
a brownian particle does not have a well-defined trajectory, and cannot posses sharply defined coordinate and coarse- 
grained velocity, as seen below. 

We saw above that the coarse-grained velocities w+j (x, t) and w_.j(x, i) are defined with respect to the ensemble 
S(x, t) of all brownian particles which at the moment t pass via coordinate vector x. On the other hand, S(x, t) is by 
itself a subensemble embedded with probability (weight) /'(x, t) into the ensemble of all realizations of the random 
coordinate vector x at the time t. Thus, both v+.j{x,t) and v^j{x,t), as well as the osmotic velocity Uj{x,t) are 
random quantities as functions of of the random configuration of N brownian particles x — {xi, ■ ■ ■ ,xn). In other 
words, this randomness enters via the context (x, t) which was chosen to define the ensemble I](x, t). 

More specifically, let us focus on the common distribution function /'(x, u; t) of 

X = (.Ti, ...,XAr) and u = (mi, mat), (27) 



N 

P(x, u; t) ^ P{^,t)Y[S{ u, - u, (x, t) ) , (28) 



where Uj (x, t) is defined by H26|) . 

One has after averaging over P(x, u; t) 



(uj) = J duj dxuj P(x, t) 6{ Uj - Uj{x, t) ) = -Tj j dxd^^P{x, t) = 0, (29) 



{{xk - {xk)) {uj - (uj) )) ^ j dxuj{x,t) P{x,t) {xk - (xfe)) = -Tj J dx{xk - {xk)) dxjP{x,t) = TjSkj- (30) 

The fact of (uj) = is natural, since uj is the (coarse-grained) velocity difference generated due to the interaction 
to the bath. It is also seen from l|3U|l that Uj of the corresponding brownian particle correlates only with its own 
coordinate. This is related to the independence of the thermal baths which act on different brownian particles. 
Applying to (|30() the standard Cauchy-Schwartz inequality one deduces: 

{{xk- (xk) ?){{u,- {u,) f)>\{{xk- {xk) ){u,- {u,) ))p = T]5k„ (31) 

which implies an uncertainty relation between the uncertainty of the coordinate and that of the osmotic velocity |29j| . 

Obviously, the regular counterpart fj{x,t) of the coarse-grained velocity does not (and should not) have any 
uncertainty relation with the coordinate. 



B. Discussion. 



Let us discuss in more detail the difference between the physical meaning of the brownian uncertainty relation 
(|31|l as compared to the standard quantum mechanical uncertainty relation ((Sj. The latter refers to two identical 
quantum ensembles £{ p{t) ) such that on the first (second) ensemble one measures x (p), and then combines into the 
single relation |SJ| the obtained dispersions. Note that the frequent interpretation of (|SJ) as referring to disturbance of 
the momentum when measuring the coordinate (and vice versa) is in general incorrect T^iS^j, since, in particular, 
the quantities (Ax'^) and (Ap^) clearly refer to different measurements done on different ensembles, rather than to a 
simultaneous measurement of the coordinate and momentum done on a single ensemble. More refined arguments and 
explicit counterexamples for this point can be found in 19, 32] . 

In contrast to this, the quantities entering into Eq. refer to consecutive measurements done on the same 
ensemble. The involved coarse-grained velocity is operationally well-defined and characterizes the change of the 
coordinate. This is contrasting to the situation with the quantum mechanical observable of the momentum which 
(in the Heisenberg representation ) in general does not fully characterize the change of the coordinate for any initial 
state; see Footnote for more details. 
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IV. BROWNIAN ENTANGLEMENT. 



Consider now two brownian particles. An analog of entanglement can be introduced in the following way. Call the 
state of two particles non-entangled (separable) if the common distribution function H28|l of the coordinates and the 
osmotic velocities can be prepared by mixing non-correlated distributions * 

P{xi,ui,x2,U2) = J dXViX)Pi{xi,ui\\)P2ix2,u2\\), 7'(A)>0, jd\V{\)^l. (32) 

We shall naturally require that the separate distributions Pi(xi,mi|A) and P2(a^2,U2|A) ensure the basic properties 
1201) of the osmotic velocity Uj . For the rest they can be completely arbitrary. The requirements (|^ 1201) make the 
basic difference as compared to the (naive) classical discussion in section FlI CI where any distribution was admissible. 
More detailed discussion on the physical meaning of H32|l . and, in particular, on its similarities and differences with 
the definition of quantum entanglement, is postponed till section IVTl 

For simplicity we choose units such that x and u have the same dimension as \/T (e.g., we measure time in units 
of square root of the damping constant; see Appendix). We take also Ti — T2 — T again for simplicity. In the same 
way as we derived (jsl I13|l . one can write the brownian uncertainty relation H31|l as 

( Am^ ) + ( Aa;^ ) > 2T, ( Am^ ) + ( Ax^ ) > 2T, (33) 

and proceed to derive 

((Aui+Au2)') + ((A2;i-A:e2)') < 4r, (34) 

as a sufficient condition for the entanglement. More general relations H14(l can also be obviously transferred to the 
brownian situation: 

( ( Aui + CAu2 f) + {{/^xi+ eAx2 f ) < 4T, (35) 
where and e can independently assume values ±1. 



A. Gaussian state of two interacting brownian particles. 



It will be useful to work out a simple model for two coupled brownian particle, where the above general concepts 
can be visualized and studied in detail. 

Consider now two harmonically interacting brownian particles with an overall potential energy 



U{xi,X2) 



gxiX2, 



(36) 



where a > and g characterizes the interaction between the particles ^ . The particles interact with independent bath 
at the same temperature T. 

Fokker-Planck equation HI8|I for the dynamics of these brownian particles can be solved directly and the outcome 
is known to be given by the following two-dimensional gaussian (provided the initial distribution was gaussian): 



VdetC 
F[xi,X2;t)= — — — exp 




0-11 (t) <7l2{t) 
0-12 (i) 0-22 (i) 



(37) 



We consider only the common distribution function of the coordinates and the osmotic velocities, and not, e.g., the common distribution 
of the coordinates and forward velocities j' (x, t), simply because due to the brownian uncertainty relation i3Ui only the former can 
lead to a non-trivial definition of entanglement. Recall in this context that the quantum uncertainty relations are necessary for the 
existence of the quantum entanglement. 
Note that the positivity of U{xi, X2) is implied by < . 
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where for simplicity we assumed that 

= (a;2W) =0, (38) 
for aU times. It is also useful to note the marginal distributions of P{xi, X2',t): 

r.2 



Pj (xj ; t) = exp 



J = 1,2. (39) 



The easiest way to obtain (Tii(t), CT22(i) and CTi2(i) is to look directly at the Langevin equations ()17|l with the 
potential H36(l : 

xi{t)^'axi{t)+gx2{t)=m{t), i2{t)+ax2{t)+gxi{t)^r]2{t), (40) 

which can be more conveniently re-written in terms of the relative coordinate r_ = {xi — a;2)/2 and the center of 
mass coordinate r+ = (xi + X2)/2: 

r±it) = -ia± g)r±{t) + ^(mit) ± 772 (t)), (41) 

and then solved directly: 

r±(t) = e-(»±9)*r±(0) + i f ds e-^''^^^ {m{t - s) ± T^2it ~ s)). (42) 

^ Jo 

Eqs. I^BII^ produce for the coarse-grained velocity difference (osmotic velocity) of the first and the second particle, 
respectively, 

T{a22Xi - (T12X2) . T (0-11X2 -gi2a:^i) 
Mi(x,t):= 2 ' U2(x,i) = 2 . (43) 

a'llCr22 — 0-12 0'lif722 — 0-12 

The sufficient condition H35|) reduces to 

( ( Aui + CAu2 )2 ) + ( ( Axi + eAx2 f ) = (^^22 + ^ii - 2(^12) ^ ^ ^ ^^^^^ < 4r. (44) 

0'licr22 — cri2 

A simpler expression is gotten for cth = 022 

r2 



fii + Co'12 



fill + ecri2 < IT. (45) 



This condition is naturally not satisfied for cri2 = 0, as the minimal value of cth -|- T^/cth over a\\ is just equal to 
2T. For the correlated situation a\2 > and we apply (|45|l for C = 1, e = — 1, while for the anti-correlated situation 
(7i2 < we use in the same way ()45|l for C = — 1, e = 1. With these we get from (|45|l the following condition 

ian~Tf <<jl2 + 2T\ai2\. (46) 

It is seen that this sufficient condition for the brownian entanglement is symmetric with respect to correlation and 
anticorrelation, and it can be satisfied if an is sufficiently close to T and \(Ti2\ is finite. Eq. (|46|l can be also satisfied 
when the particles are either correlated or anti-correlated sufficiently strongly, e.g., when |cri2| is sufficiently large and 
sufficiently close to its upper bound an. Note that (|46|) is not satisfied for T ^ 0. This is natural, since there are no 
brownian uncertainty relations in this limit. 

Thus the brownian entanglement can exist for sufficiently strongly fluctuating and/or sufficiently strongly interacting 
brownian particles. 

Let us check (|45|) more specifically with the stationary (equilibrium) state of the brownian particles which is 
established for long times provided the potential energy U{xi,X2) is (strictly) positive, and where 

Ta Tg 

"■11 = (722 = ^5 n, a\2 = 5 5-. (4('j 

~ g — 



This equality will be satisfied for two identical brownian particles if their initial conditions are the same. 
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These relations can be obtained either by directly solving H42() and then taking the limit t oo, or directly via the 
stationary distribution for the particles which is known to be Gibbsian: P{xi,X2) oc exp[-~U{xi, X2)/T]. 
Eq. (jinil reduces to 



\g\ > -1 + ^l + {a-iy. (48) 
Thus if \g\ is large enough, there can be entanglement in the equilibrium state. 

B. Entanglement for two non-interacting brownian particles. 

In the quantum case two subsystems can be entangled even if they interacted in the past, and they do not interact 
at the moment when the entanglement is tested. This is also the case with the brownian particles, as we show now. 

Return to the above system of two harmonically interacting brownian particles and assume that they did interact 
for t < but are not interacting for positive times^^. One has from 14()|l (or from 142|l ): 

=e"'"Sj(0) + ^(1-6-2'^*). <Ji2{t)=e~^'''ai2{0), J = 1,2. (49) 

We assume that the condition H45I) (with an = a22 at all times) was satisfied at the initial time, so that the 
brownian entanglement was present. For positive times and a > 0, (7i2{t) gradually disappears. This means that 
for long times the sufficient condition (|45|l for the brownian entanglement will not be valid. Clearly, no brownian 
entanglement can persists for very long times, since the particles become in this limit completely non-correlated, i.e., 
the common probaility distribution for t — > oo factorizes P{xi, X2,t) = P{xi,t)P(x2,t) (compare with our discussion 
after 

It is however to be stressed that for certain not very long times we shall surely have the condition H45(l satisfied, 
at least once it was satisfied initially. Moreover, if even this condition was not satisfied at i = 0, it can curiously get 
valid for some (not very long) times. Here is an example. For free brownian particles with a = 0, the correlation 
fi2(^) does not change in time at all, but instead the dispersion of the coordinates increases linearly with time: 
— '^ji(O) + "^Tt. Condition (|46() for the presence of entanglement now reads: 



<7niO)-T + 2Tt\ < JaU0) + 2T\ai2m- (50) 



Even if this condition was not valid at t — 0, |o'ii(0) — T\ > \/ erf 2(0) + 2r|o-i2(0)|, and additionally o-ii(O) < T, 
Eq. H50(l can still be satisfied for 

t^<t<t+, (51) 

where 



'± = 2T ■ 

Our general conclusion is that the brownian entanglement is possible if the two brownian particles do interact or 
had interacted strongly enough. The absence of interactions at present need not always to destroy the brownian 
entanglement, for some finite times it may even facilitate the suHicient conditions for its existence, as we saw above. 

V. LOCAL OSMOTIC VELOCITIES. 

In section Unl we have seen that the definition of the coarse-grained velocities H21ll22l[2Hjl is given via the ensemble 
Tilxi, X2,t) in the common context of the two brownian particles. 



Conditions 1481 are compatible with the stability requirement a > \g\, which comes from demanding positivity of U{xi, X2). Indeed, 
once \g\ satisfies 1481 . a > \g\ becomes a{2 + a) > (a — 1)^ and is satisfied provided a is not very small. 

To avoid possible misunderstanding, recall that the thermal baths acting on two particles were taken to be completely independent of 
each other, so that there is no influence via them. 
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Let us now turn to the local velocities determined with respect to the ensemble T,i{xi,t), which is obtained by 
measuring at time t only the coordinate Xi of the first particle. 

There are two completely equivalent ways for determining the coarse-grained velocities over this ensemble. The 
first way amounts to repeating definitions H21I I22|) for the ensemble J^i{xi,t) 

i^±,i{xi,t) ^ ±\im^^+o dyi— — P{yi,t ± e\xi,t), (53) 



where all the involved probability distributions contain no references on the second brownian particle. The corre- 
sponding osmotic velocity reads in complete analogy to (|26|l : 

= iy^^i{xi,t) - iy+s{xi,t). (54) 

The second way is to note that if observer 1 has no information at all from observer 2, then he effectively sums 
the ensemble T,{xi, X2, t) over all possible results of the second coordinate X2 at time t with his result xi being fixed. 
Each value of X2 during this summation is then met with probability P{x2,t\xi,t), and this results in 

l^±,l{xi,t)^ J dx2V±^l{xi,X2,t) P{x2,t\xi,t). (55) 

The equivalence between definitions H53|) and (|55|l can be established via (|24l I25|l and the Bayes formula 

Completely similar definitions can be given for the second particle, employing the ensemble T,(x2) obtained by 
measuring the second particle's coordinate only. For the corresponding osmotic component one now has from (|54l 1551 

(x, , t) = ~T,d^^ Fix, ,t), J ^1,2. (56) 
For the brownian particle in the harmonic potential, as described e.g. by H39|l . one has 

M,(a:„t)-T,^. (57) 

If now the definition 132|) of non-entanglement is applied to the common distribution function 

P{xi,fii,X2,fJ.2) = ^'l(a;i) - tll{xi))P2{x2)6{fl2 - fJ.2{x2) ), (58) 

then it is seen to be satisfied trivially: No entanglement occurs with locally defined osmotic velocities 



VI. DISCUSSION. 



Let us compare in more detail the physical meaning of the quantum mechanical entanglement versus its brownian 
analog. 

• In analogy to continuous- variable quantum entanglement, brownian entanglement is defined as a type of cor- 
relation between the coordinates and the changes of the coarse-grained velocities (osmotic velocities) of two 
brownian particles, which is impossible to reproduce by mixing non-correlated — that is, referring to each par- 
ticle separately — distributions. 

• The brownian uncertainty relation (|30|l between the coordinate and the change of the coarse-grained velocity is 
necessary for the very existence of the brownian entanglement. This is again similar to the quantum situation, 
where the analogous role of a necessary condition is being played by the quantum uncertainty relations. 



Note that for obtaining the local coarse-grained velocities i'±.i{xi,t) we had to average v±^i{xi,X2,t) over the conditional distri- 
bution P{x2,t\xi,t), and not over the unconstrained probability of the second coordinate P{x2,t). This in-equivalence can be il- 
lustrated with help of the following fact of probability theory concerning three random variables A,B,C: P(C) P{A\BC) ^ 
Ji;^ P(C|i?) P(A|B(7) = P{A\B). Would we adopt the second possibility, e.g., 11(2:1,*) = ^ dx2u{x\,X2,t) P{x2,t), this would 
lead us to an explicitly common-context (non-local) quantity, e.g., for the above example of harmonic oscillators we would get 
u{xi,t) = — — ^^^^^l^l^^^^-2 ^^-j 1 which explicitly depends on the dispersion (T22(t) of the second particle. If the two particles do not 

interact for positive times, but did interact in the past (for t < 0), u{xi, t) can still be influenced, e.g., by external fields which apply on 
the second particle for positive times (apparent or false non-locality). 
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• For checking quantum entanglement via sufficient condition H13() one needs to make measurements of non- 
commuting coordinate x^. and momentum for each quantum particle (A: = 1, 2). To this end one has to have 
two different ensembles of the particle Si and S2 (compare with our discussion in section HlFjl . In contrast, 
the brownian entanglement involves only consecutive coordinate measurements done on a single ensemble of the 
brownian pairs. 

• The momentum operator p in quantum mechanics — though being equal in Heisenberg representation to the 
time-derivative of the coordinate operator p — ni-^x — does not in general fully characterize the intuitive notion 
of "change of the coordinate for an infinitesimal time" 

In contrast, the brownian entanglement is about the coordinates and the coarse-grained velocities which do 
characterize the change of the coordinates on the coarse-grained scale of time. 

• Both quantum entanglement and its brownian analog can exist for subsystems which interacted in the past, but 
do not interact at the present. 

• In quantum mechanics the sufficient condition p3|l for entanglement involves correlations between the coordinate 
operators xi and X2, and the momentum operators pi, p2 of the subsystems Si and S2. Quantum entanglement 
can only be found after the results got via measurements on the corresponding subensembles for the subsystems 
Si and S2 are put together. The situation with the brownian entanglement is similar, since it also requires 
correlation experiments. 

• Here is finally the main conceptual difference between the quantum entanglement and its brownian analog. In 
quantum mechanics the above operators of coordinate and momentum pertain to the corresponding subsystems 
Si and S2 independently of the full system context. By this we mean that all the statistics of, e.g., pi can be 
collected via local measurements on the corresponding quantum subensemble, whether or not this subensemble 
forms a part of any larger ensemble. 

In contrast, the very definition of the coarse-grained velocities (|21II22|[2B1) involves a global (that is, depending 
on the two subsystems) ensemble I](x, i). As seen in section the purely local definition of coarse-grained 
velocities can also be given, but there will not be any entanglement on that level, for the same reason as there 
is no entanglement in other classical systems (see section III Cll . 

This conclusion on the main difference is close to the analogous conclusion of Ref . |^ , which discusses similarities 
between quantum entanglement and certain correlations in classical optics. 

VII. POSSIBILITIES OF EXPERIMENTAL REALIZATION. 

One hopes that an experimental verification of the brownian entanglement is going to be easier than that of its 
quantum analog: since Brown's discovery in 1828, various examples of brownian motion are routinely observed in many 
systems. On the other hand, the basic conditions needed for observation of the brownian entanglement amount to 
two coupled brownian particles and a resolution of the brownian motion sufficient for observing the osmotic velocities. 
Recall that on the coarse-grained time scale the osmotic velocities can be visualized as average kicks got by the 
brownian particle due to its interaction with bath particles. 

Here we discuss two experimentally studied examples of brownian motion. The pccularity of these examples is that 
the brownian motion can be detected with human eyes (without microscopes). These are thus examples of macroscopic 
brownian motion. Recall that the typical examples of brownian motion involve much smaller scales: pollen molecules 
in water observed originally by Brown had a size « 10"'^ cm, and even smaller lengths are typical for many other 
realizations of the brownian motion. 

The first experiment considers a two-dimensional circular container with elastic walls. Inside of the container 
there are motorized balls with a mass of 120 g and a diameter of 8 cm. Each ball is driven by a battery-powered motor 



This is a general point. A difference A{t) — A(0) = Jq dt^A of Heisenberg operators does not fully characterize the change in time 
of the observable A, because there are ensembles described by (time-independent in the Heisenberg representation) states \ip){ip\ for 
which {A{t) — A{0))\ili) = 0. This seems to imply that the value of A did not change at all, but this is not correct, since the above 
eigenvalue relation may be still compatible — due to non-commutativity [A(t),A(0)] ^ — with different statistics of A{t) and A(0), 
e.g., {■\p\A^ {t)\il>) ^ {xI)\A!^ {(i)\ip) . For more elaborated discussion and concrete examples see 35]. For unbound operators, such x and p, 
the above reasonings may need to be technically modified, since eigenstates of unbound operators are not normalizable. In this context 
it may suffice to require (A(t) — A(0))|i/') ~ 0. 
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inside, and moves chaotically due to elastic collisions with the walls of the container and with other motorized balls. 
The set of balls (with concentration 50 balls/m^) models the bath particles. The brownian particles are modelled with 
a set of ping-pong balls (with a mass 2 g and diameter 4 cm) which undergo random collisions with the motorized 
balls. The primary purpose of '44"! was to study polymer statistics, so that the ping-pong balls were jointed into 
a long chain by means of springs. The resulting coupled brownian motion produced a probability distribution for 
the end-to-end distance of the ping-pong chain in a close agreement with the existing theories of two-dimensional 
(self-avoiding) random walk. In this experiment Langevin equation as such was not tested directly, but there were 
several indirect reasons supporting its validity |43 |. This system would fit for observing the brownian entanglement, 
because almost all needed ingredients of this phenomenon are present. 

A curious pedagogical analog of the above experiment is the motion of a child toy called Bumble Ball • ^ plastic 
sphere about 11 cm in diameter with a symmetrical pattern of rubber knobs extending about 3 cm from its surface. 
An internal motor rotates an off-axis mass such that when activated and placed on a hard surface, a Bumble Ball 
simulates a random walk [45l |. Joining two such balls with a spring and placing them into a container creates a toy 
model which might provide at least the most rough indications of the brownian entanglement. 

The second experiment 'i^ studied a single sphere (a ping-pong ball of a mass 2.5 g and diameter 4 cm) rolling 
stochastically in an upflow of gas under flow speed 280 cm/c. The sphere rolls stochastically due to the turbulence 
it generates in the gas stream (the corresponding Reynolds number is ^40). The study focussed on the full time- 
dependent dynamics of the sphere, as captured by video imaging. It appears that the dynamics is that of the 
underdamped brownian particle in a harmonic potential. This correspondence was thoroughly checked from various 
perspectives. When looking for the brownian entanglement in such a situation, it would be necessary to go to the 
overdamped motion regime (e.g., by decreasing the mass of the sphere and by increasing the flow speed of the stream) 
and to add there a second sphere. 



VIII. CONCLUSION. 



We have uncovered the phenomenon of brownian entanglement: a correlation effect between the coordinates and the 
coarse-grained velocities of two classical brownian particles, which resembles the quantum entanglement. In contrast 
to the latter, which is presently given a fundamental status, the brownian entanglement — as the very subject of 
statistical physics — arises out of coarse-graining (incomplete description) reasons. In that respect it is similar to 
other basic relations of the statistical physics, such as the second law pj5l |. In the present situation the coarse-graining 
comes due to the time-scale separation: the evolution of the momenta of the brownian particles is very fast and 
cannot be resolved on the time-scales available to the experiment. The idea of time-scale separations is one of the 
most pertinent ones in non-equilibrium statistical physics. In a qualitative form it appears already in good textbooks 
Is^ . and has been since then formalized in various contexts and on various levels of generality 



Once there is time-scale separation between the (stochastic) motion of the coordinate and the momentum of the 
brownian particles, the operational definition of velocity via the rate of the coordinate leads to the coarse-grained 
velocity which is not equal to the real momentum. The change of the coarse-grained velocity is controlled by the 
brownian uncertainty relation. Moreover, the coarse-grained velocities appear to be contextual random quantities, 
i.e., they depend on the concrete setting of the coordinate measurement used to define them. These two aspects (un- 
certainty relations and contextuality) suffice to define brownian entanglement, similarly to quantum entanglement, as 
impossibility to prepare the common distribution of the coordinates and the coarse-grained velocities of two brownian 
particles via mixing locally independent (non-correlated) distributions referring to the two particles separately. Again 
in analogy to quantum entanglement, brownian entanglement can be witnessed via the uncertainty relations. 

In this paper we demonstrated the brownian entanglement on the simplest, exactly solvable, one-dimensional models 
of two interacting brownian particles. It should be however kept in mind that interacting brownian particles (random 
walks) are basic for several fields of modern statistical physics, such as colloids 5^ or polymers Q. One of the most 
transparent examples from the polymer science is a DNA macromolecule which consists of two interacting random 
walks (strands) As we argued, the experimental observation of the brownian entanglement might be easier than 
those of the quantum entanglement, since there are realizations of macroscopic brownian motion that are visible with 
the human eye. 



Acknowledgments 



A.E. A thanks E. Mamasakhlisov for useful discussions. The work of A.E. A is part of the research programme of 
the Stichting voor Fundamenteel Onderzoek der Materie (POM, financially supported by the Nederlandse Organisatie 



15 



voor Wetenschappelijk Onderzoek (NWO)). 



[1] N. Bohr, in Collected Works, ed. J. Kalckar (North-Holland, Amsterdam, 1985), Vol. 6, pp. 316-330, 376-377. 
[2] F. Schlogl, J. Phys. Chem. Solids, 49, 679 (1988). 

[3] L. Rosenfeld, in Ergodic Theories, Proceedings of the International School of Physics "Enrico Fermi," Course XIV, Varenna 

1960, ed. P. Caldirola (Academic Press, New York, 1961). 
[4] L.D. Landau and E.M. Lifshitz, Statistical Physics, I, (Pergamon Press Oxford, 1978). 
[5] H. Risken, The Fokker-Planck Equation, Springer- Verlag, 1984. 

[6] A.Yu. Grosberg and A.R. Khokhlov, Statistical Physics of Macromolecules, (American Institute of Physics, New- York, 
1994). 

[7] E. Nelson, Quantum Fluctuations, (Princeton University Press, Princeton, 1985). 
M. Baublitz, Jr., Prog. Theor. Phys. 80, 232 (1988); Phys. Rev. A 51, 1677 (1995). 
P. Garbaczewski and J. Vigier, Phys. Rev. A 46, 4634 (1992). 

G. Kaniadakis, Physica A 307, 172 (2002). 

[8] L. de la Peha and A.M. Cetto, The Quantum Dice: An Introduction to Stochastic Electrodynamics, (Kluwer, Dordrecht, 
1996). 

T. H. Boyer Phys. Rev. D 29, 1096 (1984); ibid., Phys. Rev. D 27, 2906 (1983); Phys. Rev. A 18, 1238 (1978). 
L.S.F. Olavo, Phys. Rev. A 61, 052109 (2000). 

R. Feynman and A. Hibbs, Quantum Mechanics and Path Integrals, (McGraw-Hill, New- York, 1965). 
L.E. Ballentine, Am. J. Phys., 54, 883 (1986). 

A. Khrennikov, J. Math. Phys., 45, 902 (2004); Md., 44, 2471 (2003); tbid., 43, 789 (2002). 

A. Peres, Quantum Theory: Concepts and Methods, (Kluwer Academic Publishers, 1993). 
D.I. Blokhintsev, The Philosophy of Quantum Mechanics, (Reidel, Dordrecht, 968). 
L.E. Ballentine, Rev. Mod. Phys. 42, 358 (1970). 
R.G. Newton, Am. J. Phys. 48, 1029 (1980). 
R. Balian, Am. J. Phys. 57, 1019 (1989). 
D. Home and M.A.B. Whitaker, Phys. Rep. 210, 223 (1992). 

B. D'Espagnat, Veiled Reality, (Addison- Wesley, MA, 1995). 

W.M. de Muynck, Foundations of Quantum Mechanics, an Empirist Approach, (Kluwer Academic Publishers, 2002). 
W.N. Itano et al, Phys. Rev. A, 47, 3554 (1993). 
R.J.C. Spreeuw, Found. Phys., 28, 361 (1998). 

A. Gatti, E. Brambilla, M. Bache, and L. A. Lugiato, Phys. Rev. Lett. 93, 093602 (2004). 

D. Collins and S. Popescu, Phys. Rev. A 65, 032321 (2 002). 

K.A. Kirkpatrick, Found. Phys. Lett., 16, 199 (2003); |quant-ph/0106072l 
O. Cohen, quant-ph/0310017 
T. Kriiger, quant-ph/0404091 

A. E. AUahverdyan, R. Bahan, and Th.M. Nieuwenhuizen, Europhys. Lett., 61, 452 (2003). 

H. F. Hofmann and S. Takeuchi, Phys. Rev. A 68, 032103 (2003). 

I. Fenyes, Z. Phys. 132, 81 (1952). 

L. de la Pena-Auerbach and A.M. Cetto, Phys. Lett. A, 39, 65 (1972). 
S. De Falco, S. De Martino and S. De Siena, Phys. Rev. Lett., 49, 181 (1982). 

S. Golin, J. Math. Phys., 26, 2781 (1985). 

R. Simon, Phys. Rev. Lett., 84, 2726 (2000); |q uant-ph/9909044| R.F. Werner and M.M. Wolf, |quant-ph/0009118l 

E. Jaynes, in The Maximum Entropy Formalism, ed. by R.D. Levine and M. Tribus, (The MIT Press, 1978). 
M. Ozawa, Phys. Lett. A, 299, 1 (2002); quant-ph/0107001 

M. Hillery, R.F. O'ConneU, M.O. Scully, E.P. Wigner, Physics Reports, 106, 121 (1984). 
Hai-Woong Lee, Physics Reports, 259, 147 (1995) 
S.-K. Ma, Statistical Mechanics, (World Scientific, Singap ore, 1985 ). 
A.E. AUahverdyan and Th.M. Nieuwenhuizen, cond-mat / 0408697| 
P. UUersma, Physica 32 27 (1966); ibid 56; 74; 90. 
A.O. Caldeira, A.J. Leggett, Ann. Phys. 149, 374, (1983). 

C. W. Gardiner, Quantum Noise, (Springer- Verlag, 1991). 
R. Landauer and J.W.F. Woo, Phys. Rev. A 6, 2204 (1972). 
L.F. Cughandolo, J. Kurchan, and L. Peliti, Phys. Rev. E 55, 3898 (1997). 
Th.M. Nieuwenhuizen, Phys. Rev. Lett. 80, 5580 (1998); Phys. Rev. E 61, 267 (2000). 
A. E. AUahverdyan and Th. M. Nieuwenhuizen, Phys. Rev. E 62, 845 (2000). 
L.S. Schulman and B. Gaveau, Foundations of Physics, 31, 713 (2001). 

W.B. Russel, D.A. Saville, and W.R. Schowalter, Colloidal Dispersions, (Cambrigde University Press, Cambridge, 1989). 
J. J. Prentis and D. R. Sisan, Phys. Rev. E 65, 031306 (2002). 
J.J. Prentis, Am. J. Phys. 68, 1073 (2000). 



16 



[45] M. Zeilik, Phys. Teach. 37, 236 (1999). 
[46] R.P. Ojha, et al, Nature 427, 522 (2004). 



APPENDIX A 



This Appendix has two closely related purposes. First, we shall apply the definition (|21l I22() of the coarse-grained 
velocities in a non-overdamped situation, that is, we shall take e in those definitions much smaller than the character- 
istic relaxation time of the (real) momenta. We shall convince ourselves that the expected answer is obtained, relating 
the coarse-grained velocities to the average (real) momentum. In addition, we shall see that the the coarse-grained 
velocity difference (osmotic velocity) (|26|l disappears in this situation, i.e., the definitions (|21II22|I become equivalent. 
Second, based on an exactly solvable situation, we shall follow in detail to the behavior of the coarse-grained velocity 
as a function of e. For simplicity we operate only with one brownian particle. 



1. 



The evolution of the common probability distribution of the coordinate x and the momentum p of the brownian 
particle is described by the following Fokker-Planck-Kramers-Klein equation [^] 

dtV = ~-d^V + dpi^p + V'{x) + jTdp)r, (Al) 

TO TO 

where to, 7 and V{x) are mass, damping constant and potential, respectively, and where 

'P{x,p,t\x',p,t'), t>t', 'P{x,p,t\x',p,t) = S{x-x')5{p-p'), (A2) 

is the conditional probability to move from [x' ,p') at time t' to {x,p) at time t. This equation corresponds to the 
Langevin equation 

mi = ~ax --1X + T]{t), {r]{t)) = 0, {T]{t)ij{t')) = 2-fT6{t - 1'), (A3) 

where —jx is the friction force, and where r]{t) is the random gaussian white noise. 
Applying the definition 1)21(1 of the coarse-grained velocity and using (|AllA2p one has 

/y — *^ 
dy ^ ' V{y,t + e\x,t) 

dydpdp' ^^V{y,p,t^ e\x,p' ,t)V{p\t\x,t) 
p 

dydp{y - x)V{p,t\x,t) —dyS{x~y) 

- / dpVip,t\x,t)^, (A4) 

where V{p, t\x, t) is the conditional probability for having the momentum equal to p at time i, provided the coordinate 
was equal to x at the same time. 
Likewise, 



X — y 

v-{x,t) =lims^+o / dy P{y,t-e\x,t) 



dydpdp —^'Pix,p,t\y,p',t- e) 

dydp{y-x) d^5{x - y) 

P(X,t) m 

dpV{p,t\x,t) — = u+{x,t). (A5) 
to 

The fact that v^{x^t) = 1^4- (a;, t) can be easily generalized to any number of interacting brownian particles. 



P(x,t) 
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Now if there are no correlations between the momentum and coordinate ^(p, i|a;,t) = 'P{p,t), then v-ix^t) = 
v^(x^t) reduce to the usual average momentum. This is the case, in particular, for the deterministic situation, where 
V{x,p,t) = 5{x-x{t))5{p-p{t)). 

Thus provided e has been taken much smaller than any relevant time-scale related to the coordinate and/or the 
momentum, the coarse-grained velocity reduces to the average momentum as it should. In this situation the osmotic 
velocity is zero. 



2. 



Now let us show on the exactly solvable situation with a harmonic potential 



ax^ 



V{x) = — , (A6) 



that if the characteristic relaxation times of the momentum 

m 



and the coordinate 



are well separated: 



and if additionally 



7 



T -1 

'X — , 

a 



(AT) 



(A8) 



Tp < r^, (A9) 



> £ > Tp, (AlO) 



one gets back the values for osmotic velocity which were obtained in the main text by means of the overdamped 
Fokker-Planck equation H18|l. 

The case with V{x) = ax^ jl can be solved either directly from |IA1|1 . or using the equivalent Langevin equation 
(|A3(I . The solution of the latter — obtained, e.g., via Laplace transformation — reads: 

x{t) = xiO)g{t) + -p{0)f{t) + - f dt'fit - t'Mt'), (All) 
TO m 

where 

m = — , m = ^ — , (A12) 



For simplicity reasons we shall put x{0) — p{0) = 0. This can be done, since these quantities are assumed to be 
independent of the noise by the very definition of the considered stochastic process. 

As the process ()A11|) is gaussian, the two-time probability distribution of the coordinate can be written down as: 

r{y, s; x,t)^^ t"^"^^ + ^22x2 + 2Ai2xy] , (A14) 



where 



d= a{t,t)a{s,s) - a^{s,t), (A15) 
/ a{s,s) CT{s,t) \ _ 1 f a{t,t) -a{s,t)\ 

[ a{s,t) a{t,t) )- d[ -<j{s,t) a{s,s) ) ' ^^^^> 
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and where <y{s,t) is the correlation function of the coordinate, 

ais, t) = {x{s)x{t)) ^-^ / dhdh f{s - h)f{t - h)S{ti - <2) - ^ / dt' fit') fit' + \s~ t\). (A17) 

Jo Jo m Jo 

Using: 

dyyViy, s\., t) ^ f y ^ ^ . . (A18) 

J Vix,t) All cr(i,i) 

one gets keeping e finite and fixed 

v+ix,t) ^fdyy^Viy,t + e\x,t)^-(^^^±^-l], (A19) 



e e \ crit,t) 

- (l - 1^ 

e \ ait,t) 



u_ix,t) ^ I dy^Piy,t~e\x,t)^^ fl-^-^], (A20) 



for any e > smaller than t: t — e > 0. 
If now e is smallest time-scale, 



V-ix, t) — v+ix, t) 



2-fT xe 



'm? ait, t) \_Jq 



dt' f\t')~fit)fit) 



+ O(e^) (A21) 



goes to zero with e ^ as was predicted above. 

As seen from (IA13|) . the overdamped limit is realized with the dimcnsionless parameter Aara/y^ being small: 

— « 1, A22 

and then the charactersitc relaxation times of the momentum and the coordinate [obtained from lAllI IA12I IA13|l ] 
are well-separated 

Tp = — w — < = ^ w — . (A23 

7 wi a UJ2 

Now if 

wii>l, wis>l, a;i|s-i|>l, (A24) 

one gets from (IA17|I 

ais, i) = ^ {e-^2\s~t\ _ ^ (A25) 

which coincides with the corresponding correlator obtained from the overdamped Fokker-Planck equation directly. 
Taking now s = t ± e in IjAlQI IA20|) and employing (|A25|) with the limit e ^ — which now should be understood 
in the context of conditions ljA23l IX24|l — one can get the expressions for v±ix, t) which can alternatively be derived 
from H24I I25|l . if we put there = 1 (only one brownian particle), fix) — —ax/^ (the linear force divided over the 
friction constant), and Pix,t) cx exp[— ^^^^ (gaussian distribution function). 



